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NOTATION 


A, B 

constants used in the Van Driest formulas for the velocity-profile transformation, see 
table 1 


C f 

local skin friction coefficient 


l 

length of plate 


M 

Mach number 


n 

fj / v \ l /n 

velocity-profile power-law exponent defined by — = 


P 

pressure 


Pr 

Prandtl number 


q 

dynamic pressure, ^ P& 2 


r 

T —T 

turbulent temperature recovery factor, _ = 0.88 

* t,e 1 e 


R x,tr 

U e 

Reynolds number based on r, , x , 

tr’ v e tr 


R 8 

U p 

Reynolds number based on momentum thickness, 6 

V 
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T 

temperature 


U 

velocity component parallel to surface 


U r 

fr 

friction velocity, a / — 

V p 


x tr 

distance from leading edge to transition point indicated by sublimation and heat transfer 
(peak heating) measurements 


y 

distance normal to flat plate 


y max 

largest height (y) in table 2 



boundary-layer thickness determined by extrapolating to U- = 1.0 the measured 

U U € 

velocity profile in power-law form, log jj- vs log y 


8 P 

boundary-layer thickness determined by linear extrapolation of pitot-pressure profiles 
to 1*0 

Me 



iii 
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boundary-layer thickness where jj~ = 0.99 


max 


boundary-layer displacement thickness, 


(' ~inr Y* 
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max 


boundary-layer momentum thickness, 


pu (^ _ uy 

o U 

Ve e 


Kb 


incremental uncertainty in quantity ( ) 
dynamic viscosity from Keyes’ formula, see table 1 
kinematic viscosity 
mass density 

scaling parameter in Baronti-Libby transformation method, see table 1 
local shear stress 


Subscripts 

adiabatic-wall conditions 
Baronti and Libby 
Coles 

boundary-layer-edge conditions 

edge of boundary-layer sublayer in Coles’ transformation method, see table 
indicated 

reservoir conditions 

total conditions (isentropic stagnation) 

wall conditions (surface) 

conditions just downstream of a normal shock wave 


Superscripts 

( ) incompressible (variable transformed to equivalent constant properties case) 

( )' reference conditions 


v 



TURBULENT BOUNDARY-LAYER VELOCITY PROFILES 


ON A N ON ADIABATIC FLAT PLATE 

AT MACH NUMBER 6.5 

Earl R. Keener and Edward J. Hopkins 
Ames Research Center 


SUMMARY 


Velocity profiles were obtained from pitot-pressure and total-temperature measurements 
within a turbulent boundary layer at one station on a large, nonadiabatic, sharp-edged flat plate at 
a local Mach number near 6.5. Local skin friction was also directly measured. Momentum- 
thickness Reynolds numbers ranged from 2590 to 8860 and wall-to-adiabatic-wall temperature ratios 
ranged from 0.3 to 0.5. Measurements were made both with and without boundary-layer trips. 

Five methods are evaluated for correlating the measured velocity profiles with the incompres- 
sible form of the law-of-the-wall and the velocity defect law. The mixing-length generalization of 
Van Driest gives the best correlation. 


INTRODUCTION 


Several investigations have been conducted at Ames Research Center to obtain measurements 
of skin-friction, heat transfer, and velocity profiles on flat plates and wind tunnel walls at hyper- 
sonic, nonadiabatic wall conditions. The purpose of the investigations has been to use the data to 
evaluate existing methods of prediction on flat plates. Correlations between the data and methods 
for predicting skin friction and heat transfer have been reported in references 1 to 5. 

This report presents an analysis of the measurements of velocity profiles in a turbulent 
boundary layer on the large flat plate model used in the investigation of turbulent skin friction, 
reported in reference 2. Local edge Mach number was about 6.5, momentum-thickness Reynolds 
number ranged from 2590 to 8860 and wall-to-adiabatic-wall temperature ratios ranged from 0.3 
to 0.5. Measurements were made both with and without boundary-layer trips. The experimental 
velocity profiles are used to evaluate five velocity-profile correlation methods. Two of the 
methods are reference-temperature methods, T-wall and T' (ref. 6); the other three methods are 
those of Coles (ref. 7), Baronti and Libby (ref. 8), and Van Driest (ref. 9). 



EXPERIMENT 


Wind Tunnel 

This investigation was conducted in air in the Ames 3.5-Foot Hypersonic Wind Tunnel, in 
which cold air was passed through an alumina storage heater system and heated to stagnation 
temperatures ranging from 684° to 1089° K. The Mach 7.4 contoured nozzle was used. Stagnation 
pressure was varied from 39 to 66 atm and the testing time was between 2 and 3 minutes depend- 
ing on the pressure level. 


Model and Instrumentation 

The model was a sharp-edged flat plate, 122.0 cm long, 45.7 cm wide, and 2.5 cm thick, 
firmly supported by a massive undercarriage instrument housing rigidly attached to the strut of the 
sting support system (fig. 1 (a)). A dimensional sketch of the flat plate is presented in figure 1 (b). 
The boundary layer was surveyed from above the plate at a station 99.6 cm behind the leading 
edge using a movable survey complex having pitot-pressure, static-pressure, and total-temperature 
probes (fig. 2). The total-temperature probe was a thermocouple surrounded by a single shield 
vented at the top behind the thermocouple. A single shield was used because it could be made 
small and, thereby, minimize possible interference effects from interaction of the probe shock and 
the boundary layer. Skin friction was measured by a floating-element balance (fig. 3). (See ref. 2 
for additional details.) Surface pressures and temperatures were measured along the model center- 
line. Wind-tunnel total pressure and total temperature were measured in the tunnel reservoir. 


Test Conditions 

The sharp-edged flat plate was mounted with its test surface 3° to the windward. (This angle 
was required for another investigation and could not be varied.) The resulting local edge Mach 
number was about 6.5 for a free-stream Mach number of 7.4. Local unit Reynolds number was 
varied from about 5 to 1 3 million per meter by varying tunnel reservoir pressure. Wall-to-adiabatic- 
wall temperature ratio was varied from 0.3 to 0.5 by varying tunnel reservoir temperature. Nearly 
isothermal wall conditions (within 30° K) existed along the centerline of the plate as verified by 
thermocouple measurements. The increase in surface temperature during a run was less than 5 
percent of the total temperature. Measurements were made with and without boundary-layer 
trips. 


DATA REDUCTION 


Local flow conditions at the edge of the boundary layer on the flat plate (99.6 cm behind the 
leading edge) were calculated from the measured surface static pressure (p w ), boundary-layer- 
edge pitot pressure (p t 2 ) e , and the reservoir total temperature (7^ 0 ). The compressible flow 
relations in reference 10 were used in the calculations, which included 1 corrections for calorically 
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imperfect gas effects as functions of total temperature. Keyes’ formula (ref. 1 1) was used for vis- 
cosity. The boundary-layer velocity profiles were calculated from the measured pitot and surface- 
static pressures (assuming that p = p w ) and from an assumed linear variation of total temperature 
with velocity (i.e., U/U e = ( T f - T w )j{T t e - T w )), which is identical to the well-known Crocco 
temperature distribution for Prandtl number of unity. Calculated, rather than measured, total 
temperatures were used because measured 7^ was not available for all of the data; however, 
where available, the measured T t agreed closely with the calculated temperature distribution 
(as shown in the discussion section). Additional details of the data reduction procedures are 
given in reference 2. 

Six velocity profiles were selected from the available test conditions (ref. 2) using the crite- 
ria of turbulent flow as indicated by the skin friction results, and constant p t Q within 2 N/cm 2 
during each boundary-layer survey. Calculated results are tabulated in table 2. 


ACCURACY 


The estimated probable uncertainties of the recorded and calculated quantities were deter- 
mined in reference 2 and are given below: 


M 

e 

±3 percent 

(p l.2>e- V T w 

±2 percent 

Pe , T t e , 8*, e, t w 

±5 percent 

R 9 

±8 percent 

y 

±0.013 cm 

u 

±2 percent 


U 

e 


VELOCITY-PROFILE CORRELATION METHODS 


Five correlation methods were selected for evaluation. Each method contains functions for 
correlating compressible velocity profiles with a generalized incompressible profile. Correlation 
functions for each method are presented in table 1 on the basis of the law-of-the-wall, 


JJ 

U 

T 
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and the velocity defect law 


U-U 
e_ 

U 



Two of the methods (T-wall and T ' (ref. 6)) follow a reference-temperature approach that was 
originally developed for the correlation of compressible skin friction. Two methods (Coles, ref. 7, 
and Baronti and Libby, ref. 8) follow a corresponding-stations approach developed for the trans- 
formation of the complete compressible turbulent boundary layer to a corresponding incompres- 
sible boundary layer. The last method (Van Driest (ref. 9)) uses the mixing-length approach 
developed for the correlation of compressible skin friction. (Use of either the Prandtl mixing 
length (ref. 9) or the Von Ka'rman mixing length (ref. 12) results in identical functions for 
eqs. (25) and (29) of table 1.) 


RESULTS AND DISCUSSION 


Before presenting the analysis of the velocity profiles, several other measurements pertinent 
to the analysis are discussed. These are skin friction, transition location, surface pressures, and the 
total-temperature and pitot-pressure surveys. 


Local Skin Friction 

Local skin friction (surface shear stress) appears as a parameter in all of the correlation for- 
mulas for the velocity profiles (table 1). Since measured (rather than theoretical) skin friction is 
used in this analysis, the measured skin-friction coefficients from reference 2 are again presented 
in figure 4 as a function of the measured momentum-thickness Reynolds number. The solid sym- 
bols correspond to the six boundary- layer profiles in table 2. The cross-hatched bands represent 
the skin friction by the theory of Van Driest (II) 1 (ref. 12). The band width represents the effect 
of the experimental range of T w /T aw . The dashed curve represents the expected level of the data 
for laminar flow according to the theory of Rubesin and Johnson (ref. 13). 

It was a conclusion of reference 2 that the theories of Van Driest (II) (ref. 12) and of Coles 
(ref. 7— not shown in fig. 4) give the best predictions of skin friction for turbulent flow over a flat 
plate for these test conditions. 


Boundary-Layer Transition 

The sublimation technique was used to determine the distance from the leading edge of the plate 
to the end of transition to verify that turbulent flow existed at the measuring station. The transi- 
tion Reynolds numbers for natural transition found by using this technique are presented in 

1 The (II) refers to the second skin-friction theory of Van Driest (1954, ref. 12), which uses the Von Karman mixing length in 
the Prandtl shear-stress; the Van Driest (I) theory (ref. 9) uses the Prandtl mixing length. 
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figure 5 as a function of unit Reynolds number. The sublimation results show that R x tr appears 
to increase with unit Reynolds number (so-called “unit Reynolds number” effect) much the same 
as the prediction for the end of transition by the semiempirical method of Deem and Murphy (refs. 
14 and 15) for 0.013 cm leading-edge thickness. Figure 5 also shows that the sublimation results 
generally correlate with transition Reynolds numbers determined from peak heating measure- 
ments obtained on a heat-transfer flat-plate model of the same size tested in the same facility 
(data of T. E. Polek, unpublished). 


Surface Pressure Distributions 

Figure 6 presents the measured surface pressures along the centerline of the model for the 
conditions of boundary-layer trips off and trips on. Location of the survey station (x/l = 0.817) 
is shown by an arrow. The pressure ranged from 0.20 to as low as 0.063 N/cm 2 . (The estimated 
probable error, ep w — ±0.0048 N/cm 2 is indicated on the figure for comparison.) The pressures 
are not as uniform as expected for a flat plate and the largest variation in surface pressure occurred 
for the highest surface pressure with boundary-layer trips on. Calculated inviscid pressures for a 
two-dimensional flow-deflection angle of 3° and Moo~ 2.37 (solid points plotted at the survey 
station) are generally within experimental accuracy of the measured pressures near the survey 
station. 


Boundary-Layer Total-Temperature Profiles 

The total-temperature probe (fig. 2) was calibrated from measurements of the probe outside 
of the boundary layer where the correct total temperature was taken to be that measured in the 
reservoir (T t Q ). The calibration followed a method similar to that used by Winkler (ref. 16) in 
which the temperature recovery factor (T t ^ — T e )j{T t — T e ) was shown to correlate with a param- 
eter based on the Nusselt number of the stagnated flow inside the probe and to be independent of 
Mach number. Winkler’s correlating parameter was rewritten as (p t 2 )(7^)~ 7 / 4 and in figure 7 the 
calibration data are shown as a function of this parameter. Measurements were obtained over a suf- 
ficient range of tunnel total pressure to cover the range of flow conditions within the boundary- 
layer survey, as shown by the arrows. The calibration data at (p t 2 )(7^)~ 7 ^ 4 > 6X10~ 6 were 
different for the two free-stream total-temperature levels of 780° and 1 050 K. 

Boundary-layer total-temperature profiles measured for three of the six velocity profiles pre- 
sented (table 2(a), (b) and (d))are shown in figure 8. Static temperatures calculated from T t and 
M are presented as T/T e . The measured temperatures agree well with the calculated Crocco tem- 
peratures for each data point (using the linear relationship, written as a function of T w and M, 
see ref. 2). The linear relationship is further verified in figure 9 for the conditions shown for a non- 
adiabatic flat plate. 


Boundary-Layer Pitot-Pressure Profiles 

The six pitot-pressure profiles from table 2 are presented in figure 10. Three profiles are for 
natural transition and three are for forced transition using boundary-layer trips. Results show that 
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the trips cause thickening of the boundary layer. The figure also shows one method of determining 
a turbulent boundary-layer thickness, obtained by linear extrapolation of the pitot-pressure pro- 
files. The resulting values of 8 p> shown by the arrows, are recorded in table 2. 


Boundary-Layer Velocity Profiles 

The boundary-layer velocity profiles are presented in figure 1 1 . There are no obvious probe- 
interference effects near the wall and all profiles are typical of well developed turbulent profiles. 
The boundary-layer thickness (8 U ), where the value of UlU e = 0.99, is shown by the tick marks and 
recorded in table 2. 

Figure 12 presents a log-log plot of y vs ( U/U e ) for the velocity profiles of figure 1 1 to 
test the power-law relationship. 


§■ =/0) 1/n 

u e . 


The data follow the linear power-law fairings quite closely from U/U e =1.0 to the data point 
closest to the wall. The boundary-layer thickness (8 }l ), determined from the intersection of the 
power-law fairing with UlU e = 1 .0, is shown by the arrows and the values are recorded in table 2. 

In figure 1 3 the velocity-profile power-law exponent is presented as a function of Reynolds 
number Rq. Included is an adiabatic-wall curve (Fenter, ref. 17) that approximately represents 
existing experimental results for subsonic and supersonic Mach number. (Fenter’s analysis indicated 
that the power-law factor is not a function of Mach number.) These results clearly do not follow 
the trend with Reynolds number of Fenter’s curve. Instead, the value of n is about double that of 
the subsonic/supersonic value of about 5 at Rq — 2500 and approaches the subsonic/ supersonic 
value of about 7 atRg — 10 4 . Recently, Johnson and Bushnell (ref. 18) summarized a large number 
of existing measurements, including the present results, and reported that the power-law expo- 
nents at M > 4 are often as high as 10 or more (termed overshoot) in the low Reynolds number 
range for flat plates, cones, and hollow cylinders but that data from wind-tunnel walls often do not 
show this effect. 

Results from four other investigations are also shown in figure 13. Boundary-layer measure- 
ments on the walls of the Ames 3.5-Foot Hypersonic Wind Tunnel at M e = 7.4 (ref. 19) and the 
Ames 8- by 7-Foot Supersonic Wind Tunnel (unpublished) are represented by two identified lines, 
which lie within 10 percent of Fenter’s empirical curve. (Skin-friction measurements from the wall 
tests were reported in refs. 2, 19, and 20). Results from the Naval Ordinance Laboratory boundary- 
layer channel at M = 5 (ref. 2 1 ) are represented by another line, which also lies within 1 0 percent 
of the Fenter curve, although at a slightly different slope. A wind-tunnel wall measurement by Hill 
(ref. 22) at M = 9 and low Reynolds number (Rq = 2000) lies on Fenter’s curve. Finally, the 
cross-hatched area represents a more recent investigation (ref. 3) extending the results to 
Rq — 1.8X1 0 4 using a flat plate that was very similar to that used in the present investigation. 
The power-law exponent appears to follow close to the subsonic/supersonic variation for 
^>7X105. 
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Boundary-Layer Thickness 

The edge of a turbulent boundary layer is difficult to determine accurately from experimental 
profiles. Neither 5* nor 0 require an accurate value for 5. However, the velocity-defect cor- 
relations require a more accurate determination of boundary-layer thickness. Three methods for 
determining 5 have been described in the discussion of the pitot-pressure profiles (8^), the 
velocity profiles (5 W ), and the power-law exponent (8 n ). Tabulated results from table 2 show 
that 8 U is close to 8 n whereas 5^ is usually much larger. The power-law method was chosen for 
use in the equations for y/8 in table 1 because this method gives consistent and rational values 
and because the power law is a characteristic expression for the outer part of the velocity profile. 


Velocity-Profile Correlation Methods 

The five velocity-profile correlation methods given in table 1 were selected for evaluation 
using the velocity profiles listed in table 2. The transformed velocity profiles are compared in 
figures 14 to 19 on the basis of the law-of-the-wall and the velocity-defect law, which have been 
shown to correlate incompressible profiles. Each method is evaluated by observing how well the 
compressible (variable-density) velocity profiles are transformed to an incompressible (constant- 
density) profile. Coles’ incompressible listing (table 3), determined empirically (ref. 23), is used 
in this analysis to represent the incompressible profiles. 

The deviation of the transformed profiles from Coles’ curve for the law-of-the-wall can be 
expressed in terms of an equivalent percentage deviation in skin friction (j w ), which appears in 
the correlating parameters. The data were transformed using measured r w rather than the theo- 
retical t w that would be obtained from the comparable skin-friction law for Cr = f (Rq ) for 
each method (ref. 2). On the other hand, a value of t w could have been chosen that would have 
forced the data to agree with Coles’ incompressible law-of-the-wall curve, and the resulting t w 

could have been compared to the measured value. The latter method, known as theClauser tech- 
nique for obtaining skin friction from velocity profiles, is described in reference 24. For discus- 
sions of previous analyses of incompressible and compressible velocity profile results, see ref- 
erences 25 and 26. 

Figures 14 to 18 show the results of using the five correlation methods on three of the six 
measured velocity profiles, one at each temperature ratio tested ( T w /T aw = 0.32, 0.43, and 0.5 1 ), 
obtained with boundary-layer trips on. 

T-wall method — The simplest correlation method is the wall reference-temperature method. 
This method is of interest because it has been used successfully to correlate the effect of compres- 
sibility on laminar skin friction and velocity profiles for adiabatic flat plates at Mach numbers up 
to 5. It is readily apparent in figure 1 4 that the T-wall method does not correlate the data since the 
three profiles do not collapse to a single curve in the law-of-the-wall transformation. The velocity- 
defect transformation does not correlate any of the profiles with the incompressible curve. 

T' method— Another simple correlation method is to evaluate the local fluid properties at a 
/’'reference temperature, similar to the reference-temperature correlations for skin friction. The T' 
method was first shown by Rubesin and Johnson (ref. 13) to extend the correlation of laminar 
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skin friction to Mach numbers much higher than 5. The reference temperature expression (eq. (6)) 
of Sommer and Short (ref. 6) was chosen to represent this type of transformation. Figure 15 shows 
that the T' method correlates the law-of-the-wall profiles better than the T-wall method. Expressed 
in terms of skin friction (r w ), the transformed profiles are within 8 percent of Coles’ incompres- 
sible curve. For example, an 8 percent higher r w than that measured would put the data for 
T w lT m = 0.32 and 0.43 onto Coles’ curve. Again, none of the velocity-defect profiles are cor- 
related with the incompressible curve. 

Method of Coles (ref. 7) — Figure 16 presents the results of the transformation of Coles’, 
developed from his law of corresponding stations. The transformed profiles lie parallel to the 
incompressible law-of-the-wall curve but are displaced as much as 1 3 percent in t w . It is interesting 
to note, however, that Coles’ skin-friction equation (ref. 2) predicts the measured skin friction 
within 5 percent Fvidentlv. the skin friction predicted by the transformation functions for skin 
friction can be different from the skin friction obtained from the velocity profiles using the func- 
tions derived from the same method. Correlation of the velocity-defect profiles with the incompres- 
sible curve is poor. 

Method of Baronti and Libby (ref 8) — In figure 17 the results are presented for the trans- 
formation of Baronti and Libby, which is an extension of Coles’ method. The method correlates 
the effect of temperature by collapsing the law-of-the-wall profiles; however, the transformed 
profiles have a different slope than the incompressible curve. These transformed profiles are 
similar to those shown by Baronti and Libby at high supersonic Mach numbers and by Bertram 
et al. (ref. 25) and at hypersonic Mach numbers by Watson and Cary (ref. 26). Skin friction can be 
obtained from the velocity profiles within 10 percent of the measured values if an abscissa 
(U T y/v) of about 200 is used to correlate with Coles’ curve. Correlation of the velocity-defect 
profiles with the incompressible curve is very poor; a similar result with this method was reported 
in reference 8. 

Method of Van Driest (ref 9) — The final correlation method considered (fig. 18) is that of 
Van Driest, developed from mixing-length theory. The method gives the best correlation of the 
law-of-the-wall profiles (within 5 percent of Coles’ incompressible curve). In addition, it is the only 
method that transforms the velocity-defect profiles close to the incompressible curve. Similar 
results are obtained (fig. 19) with the Van Driest method using the three velocity profiles in table 2 
with boundary-layer trips off (natural transition). The velocity-defect transformation shows an 
effect of Reynolds number in both figures 18 and 19. The correlation deteriorates with decreasing 
Reynolds number; however, this agrees with the power-law-exponent overshoot at the lower 
Reynolds numbers (fig. 13). The correlations of the Van Driest law-of-the-wall transformations at 
all Reynolds numbers presented are compatible with the correlation of skin friction (fig. 4) for the 
Van Driest method. The correlations indicate that skin friction can also be determined from exper- 
imental velocity profiles (Clauser technique) within ±5 percent. 

The results of figures 1 8 and 1 9 confirm that the mixing-length concept can be used to analyze 
some of the characteristics of compressible turbulent boundary layers at Mach numbers up to at 
least 7. For example, Maise and McDonald (ref. 27) chose the Van Driest generalized velocity dis- 
tribution to calculate the effect of compressibility on the shear-stress, mixing-length, and eddy- 
viscosity distributions at Mach numbers up to 5. Mathews, Childs, and Paynter (ref. 28) used the 
Van Driest generalization in Coles’ universal wall-wake equation to calculate velocity profiles for 
cases with pressure gradient and shock-wave boundary-layer interaction. 
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CONCLUDING REMARKS 


Velocity profiles were obtained from pitot-pressure and total-temperature measurements 
within a turbulent boundary layer at one station on a large (1.22 m length) nonadiabatic, sharp- 
edged flat plate at a local Mach number near 6.5. Local skin friction was also directly measured. 
Momentum thickness Reynolds numbers ranged from 2590 to 8860' and wall-to-adiabatic-wall 
temperature ratios ranged from 0.3 to 0.5. Measurements were made both with and without 
boundary-layer trips. 

The sublimation technique was used to determine the location of transition. The results gen- 
erally correlate with the location of peak heating. Measured total temperatures in the boundary 
layer follow the Crocco linear variation with velocity. Boundary-layer trips produce no obvious 
distortion effects on the velocity profiles. The velocity-profile power-law exponent does not follow 
the normal subsonic/supersonic trend of increasing with increasing Reynolds number. Instead, the 
exponent is about double the subsonic/supersonic value of about 5 at low Reynolds numbers near 
peak skin friction (Rq — 2500) and approaches the subsonic/supersonic value of about 7 at a 
momentum-thickness Reynolds number of about 10 4 . 

Five methods are evaluated for correlating the measured velocity profiles with the incompres- 
sible form of the law-of-the-wall and the velocity defect law. The Van Driest method gives the best 
correlation for both the law-of-the-wall and the velocity-defect law. Law-of-the-wall transformations 
are satisfactory at all Reynolds numbers presented, which agrees with the good correlation of skin 
friction reported previously for the Van Driest method (ref. 2). This result indicates that skin fric- 
tion can be determined from experimental velocity profiles (Clauser technique) within ±5 percent, 
using the Van Driest law-of-the-wall transformation. The velocity-defect profiles correlate approxi- 
mately with the incompressible profile at the highest momentum-thickness Reynolds numbers near 
8000. The correlation deteriorates with decreasing Reynolds number, a result that agrees with the 
power-law-exponent results. 


Ames Research Center 

National Aeronautics and Space Administration 

Moffett Field, California, 94035, March 23, 1972 
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TABLE 1 TRANSFORMATION FUNCTIONS FOR VELOCITY PROFILES 


J7 

O 


UfU e 


Wall reference temperature 

T reference temperature 

Coles (ref 7) 

Baronti and Libby (ref. 8) 




u/u e 








rz 

(10) 







y c r.c 12 







where 

c 

f ' C Ww ap e * 

(in 

where 

viu e 

V 2 

(18) 




(12) 


c .'*JLc 

fM P w d w ofi e V 

(19) 


Van Driest (ref. 9) 


uiu e 

where 

T’tT = | + 0.035 All 


•“Of-) 


"■ayym 

v T th w '-W)& 


l(rcf 6) (6) 


,y 0M, efe c ) 


U/U given as a function of U y/v in 
table 3 (taken from ref. 23) r 

equation (14) is based on aCrocco tem- 
perature distribution for unit Prandtl 
number. 

is obtained from equations ( J I ) to 
< 1 4) by iteration. 


°Kr;[ 

c jm uaa ! J (20 ) 
r T ‘ [(v) * (' + 0M ‘ "»•« 

- (I0.6)’j (21) 


C, is obtained from eqs. ( 1 9) through 
(21 ) by iteration. 


♦*" (r^=) 

YB 2 +4/t 2 / 


sin " 1 

Z ^(U/Uj-B 


. vV+4/1 2 . 




U y U a ( (C, T \ 

-f iy-y f)' 

(2) 

H IH)- 

(7) 


(15) 

U f o „ 

2 ?/„ 

(22) it | 

IV 

(A 

if) 

(28) 








i 

sin -1 

2/1 2 (tf/I/J 


v-&. (wkJ-i 

<3) 

(wu c ) -1 

(8) 

(u/U p )-l 


W - 1 

(23) •* 4/ /J )(V r c ) 


'V + 4/1 2 


u - fWV 

: WWy 

j C f.c n 

(16) 

y Cf'HL 11 




(29) 








U’ + 4„’) 



f r 

(4) 

y 

(9) 

/>.)* 

Jo 

(17) 

f ( p/p e) dy 
J o 

(24) 

3L 


(30) 



6 n 

j o ( PlPc ) dV 





NOTE: Experimental values for (y were used in above functions rather than theoretical values 
Keyes’ viscosity formula (ref. 1 1 ) used incite calculations is 

V= 1.49X10“® /r(l +-^X 10 r ) (31) 


T = temperature, * K 



TABLE 2.- BOUNDARY-LAYER SURVEYS 


M = 6.21 

e 


T. = 1 028° K 

t,e 


T = 324° K 
w 

T = 123° K 

e 

T /T =0.34 
w aw 


5 m = 1.02 cm 


5 = 1.30 cm 

P 

8 = 1 . 1 2 cm 

n 

n = 10.9 

Data symbol o 


p (ao = 394 N/cm 2 
( P ,^ = 5.8°N/cm 2 
P e = 0.1 16 N/cm 2 

(P e U J V e ) X 1 0_6 = 5 - 1 8 m_1 

U e = 1383 m/sec 
5* = 0.559 cm 
0 = 0.0498 cm 
H= 11.2 


= 2,590 
Cy = 0.00156 
B.L. trips off 


cm 

M/M e 

777, 

u/u e 

T t !T te 

measured 

0.24 

0.5651 

2.225 

0.8429 

0.907 

.27 

.5757 

2.179 

.8498 

.913 

.39 

.6316 

1.950 

.8821 

.929 

.51 

.6865 

1.753 

.9090 

.948 

.68 

.7682 

1.501 

.9413 

.963 

.74 

.7882 

1.447 

.9480 

.968 

1.08 

.9245 

1.135 

.9848 

.987 

1.26 

.9630 

1.063 

.9929 

.987 

1.28 

.9647 

1.060 

.9933 

.986 

1.41 

.9807 

1.032 

.9964 

.989 

1.68 

.9937 

1.010 

.9988 

.995 

1.85 

.9966 

1.006 

.9994 

.996 

2.42 

1.0000 

1.000 

1.0000 

1.000 


I 



TABLE 2.- BOUNDARY-LAYER SURVEYS - Continued 


(b) 


M =6.39 

e 

7\ = 1089° K 
t,e 

T = 326° K 
w 

T = 125° K 

e 

T /T =0.32 
w aw 


5 W = 1.37 cm 


p t = 666 N/cm 2 

(^, 2 ) =9 - 87N / cm2 

=0.186 N/cm 2 

(p e U/g X 10-6= 8.43 m-1 

U e = 1432 m/sec 

5* = 0.698 cm 


5^ = 1 .60 cm 

5 = 1.52 cm 

n 


6 = 0.0625 cm 
H= 11.1 


n = 7.5 


R q = 5,280 


Data symbol □ 


C f = 0.00122 


B.L. trips on 


y, 

cm 

M/M e 

1 

u/u e 

"V r „ 

measured 

0.17 

0.4509 

2.819 

0.7565 

0.854 

.28 

.5071 

2.524 

.8054 

.873 

.40 

.5480 

2.326 

.8356 

.889 

.47 

.5918 

2.128 

.8634 

.902 

.52 

.6122 

2.043 

.8750 

.913 

.64 

.6620 

1.851 

.9006 

.927 

.68 

.6792 

1.789 

.9085 

.931 

.77 

.7147 

1.670 

.9235 

.940 

.91 

.7626 

1.523 

.9412 

.955 

.94 

.7789 

1.477 

.9466 

.958 

1.01 

.8005 

1.419 

.9534 

.963 

1.06 

.8234 

1.360 

.9602 

.970 

1.12 

.8478 

1.301 

.9669 

.973 

1.25 

.8990 

1.187 

.9796 

.985 

1.30 

.9128 

1.159 

.9827 

.989 

1.46 

.9550 

1.078 

.9916 

.998 

1.78 

1.0000 

1.000 

1.0000 

1.000 
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TABLE 2.- BOUNDARY-LAYER SURVEYS - Continued 


(c) 


M 


t,e 


w 


T IT 
w' aw 


= 6.42 
= 764° K 

= 3 1 2° K 
= 85° K 
= 0.45 
= 1.22 cm 
= 1.50 cm 
= 1 .24 cm 


n = 9.4 

Data symbol 0 


P t oo = 385 N/cm 2 
( p t , 2 ) e = 6-09 N/cm 2 
„ =0.1 14 N/cm 

(^A) xl0 ~ 6=912m_I 

U =1185 m/sec 

e 

5* = 0.658 cm 
d = 0.05 18 cm 
H= 12.7 
R Q = 4,720 
Cf= 0.00123 
B.L. trips off 


y, 

cm 

M/M e 

T/T e 

u/u e 

0.21 

0.4872 

2.749 

0.8079 

.34 

.5471 

2.414 

.8501 

.48 

.6073 

2.121 

.8844 

.59 

.6600 

1.899 

.9094 

.72 

.7192 

1.682 

.9326 

.82 

.7640 

1.538 

.9475 

.97 

.8198 

1.381 

.9633 

1.09 

.8745 

1.246 

.9764 

1.21 

.9114 

1.166 

.9842 

1.28 

.9334 

1.122 

.9885 

1.48 

.9743 

1.045 

.9958 

1.61 

.9873 

1.022 

.9980 

1.72 

.9936 

1.011 

.9990 

1.81 

.9983 

1.003 

.9997 

1.99 

1.0000 

1.000 

1.0000 



TABLE 2 - BOUNDARY-LAYER SURVEYS - Continued 



n = 232 N/cm 2 

tpo 

(P,, 2 ) c = 3 - 7 ON/c- 2 

= 0.069 N/cm 2 
(pU e /v e )Xl0- 6 =5.34m-i 

U ' = 1 200 m/sec 
5* = 0.919 
9 = 0.0744 
H= 12.4 
R 0 = 3,960 
cy = 0.00125 


y, 

cm 

M/M e 

T/T e 

0.21 

0.4541 

2.937 

.34 

.4954 

2.688 

.48 

.5418 

2.433 

.58 

.5809 

2.237 

.71 

.6245 

2.040 

.82 

.6594 

1.897 

.97 

.7043 

1.730 

1.17 

.771 1 

1.515 

1.23 

.7836 

1.478 

1.48 

.8513 

1.300 

1.72 

.9135 

1.161 

3.82 

.9369 

1.114 

2.11 

.9770 

1.040 

2.21 

.9848 

1.026 

2.67 

1.0000 

1.000 


C//U e 

0.7782 

.8123 

.8451 

.8688 

.8918 

.9081 

.9264 

.9491 

.9528 

.9708 

.9845 

.9890 

.9962 

.9975 

1.0000 


B.L. trips on 

T t /T te 

measured 

0.876 

.894 

.914 

.928 

.944 

.956 

.973 

.990 

.994 

1.006 

1.006 

1.010 

1.006 

1.005 

1.000 
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TABLE 2.- BOUNDARY-LAYER SURVEYS - Continued 


(e) M = 6.50 
e 

T t = 684° K 

t,e 

T w = 318” K 

T = 73° K 

e 


T IT 
w aw 


0.51 


5 W = 1.35 cm 

6 = 1.52 cm 

P 

5 = 1.35 cm 

n 

n - 8.7 


Data symbol A 


p = 452 N/cm 2 
L°° 

( P L 2 ) e =7 - 31N/cm2 

P g =0.133 N/cm 2 

( p e*V I 'e ) X1Cr6 = 13 - 4m_1 

(/ =1119 m/sec 

e 

8* = 0.704 
0 =0.0518 
//= 13.6 
= 6,910 
Cy= 0.00106 
B.L. trips off 


y, 

cm 

M/M e 

T/T e 

u/u e 

0.14 

0.4208 

3.283 

0.7624 

.28 

.4822 

2.857 

.8149 

.40 

.5282 

2.573 

.8472 

.55 

.6069 

2.157 

.8914 

.64 

.6462 

1.980 

.9093 

.76 

.7012 

1.762 

.9308 

.89 

.7617 

1.556 

.9502 

1.02 

.8179 

1.393 

.9652 

1.16 

.8827 

1.232 

.9797 

1.25 

.9179 

1.155 

.9865 

1.41 

.9559 

1.079 

.9931 

1.55 

.9777 

1.039 

.9966 

1.63 

.9892 

1.019 

.9984 

1.68 

.9917 

1.014 

.9988 

1.95 

1.0000 

1.000 

1.0000 




TABLE 2.- BOUNDARY-LAYER SURVEYS - Concluded 


(f) M = 6.50 

e 

p. = 452 N/cm 2 

t,oo 

T = 689° K 

t,e 

(P t £ =7.29 N/cm 2 

T w = 3 1 8° K 

P e = 0.133 N/cm : 

T = 74° K 

e 

(W e ) xl0r ‘- 131 

T /T =0.51 
w aw 

U e - 1 1 23 m/sec 

8 U = 1.65 cm 

6* = 0.909 

6 =1.85 cm 

P 

0 = 0.0676 

8 = 1.52 cm 

n 

7/= 13.5 

n = 6.7 

Rg = 8,860 

Data symbol C\ 

^=0.00100 


B.L. trips on 


y, 

cm 

M/M e 

T/T e 

u/u e 

0.12 

0.3696 

3.662 

0.7072 

.25 

.4190 

3.287 

.7598 

.36 

.4486 

3.077 

.7867 

.51 

.5140 

2.653 

.8371 

.60 

.5472 

2.461 

.8584 

.74 

.5935 

2.220 

.8841 

.86 

.6405 

2.003 

.9064 

.98 

.6922 

1.794 

.9271 

1.13 

.7483 

1.598 

.9460 

1.26 

.8054 

1.426 

.9619 

1.37 

.8457 

1.320 

.9716 

1.49 

.8895 

1.216 

.9810 

1.76 

.9646 

1.063 

.9945 

2.25 

.9946 

1.009 

.9992 

2.52 

1 .0000 

1.000 

1 .0000 




TABLE 3.- INCOMPRESSIBLE LAW OF THE WALL AND VELOCITY DEFECT LAW VALUES 

[from ref. 23] 


Law of the wall 

Velocity defect law 

U T J 

u 

y 

u-u 

e 

V 

u 

T 

8 

w 

r 

0 

0 

0.010 

-14.31 

1 

.99 

.015 

-13.30 

2 

1.96 

.020 

-12.58 

3 

2.90 

.025 

-12.02 

4 

3.80 

.030 

-11.57 

5 

4.65 

.040 

-10.85 

6 

5.45 

.050 

-10.29 

7 

6.19 

.060 

-9.83 

8 

6.87 

.080 

-9.11 

9 

7.49 

.100 

-8.56 

10 

8.05 

.150 

-7.54 

12 

9.00 

.200 

-6.70 

14 

9.76 

.250 

-6.00 

16 

10.40 

.300 

-5.37 

18 

10.97 

.350 

-4.79 

20 

11.49 

.400 

-4.25 

24 

12.34 

.450 

-3.73 

28 

12.99 

.500 

-3.23 

32 

13.48 

.550 

-2.76 

36 

13.88 

.600 

-2.31 

40 

14.22 

.650 

-1.89 

44 

14.51 

.700 

-1.50 

50 

14.87 

.750 

-1.14 

60 

15.33 

.800 

-.82 

80 

16.04 

.850 

-.53 

100 

16.60 

.900 

-.29 

150 

17.61 

.950 

-.10 

200 

18.33 

1.000 

0 

300 

19.34 



400 

20.06 



500 

20.62 



600 

21.08 



800 

21.79 



1000 

22.35 



1500 

23.36 



2000 

24.08 



3000 

; 25.09 



4000 

25.81 



5000 

26.37 



6000 

26.83 



8000 

27.54 



10000 

28.10 
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@ Skin - friction balance 

(§) Boundary - layer survey complex 

(C) Preston tubes (ref. 5) 

( 0 ) Boundary - layer trips 

(a) Mounted in the Ames 3.5-Foot Hypersonic Wind Tunnel. 


Figure 1 Sharp-edged flat plate. 






(b) Geometry. 


Figure 1 Concluded. 



1.60 


L.E. of pitof probe 







Solid points: Six selected velocity profiles, table 2 
Flagged points : B.L. trips on 


Turbulent B.L. theory 
(Van Driest H , ref. 12) 



\ ( T w / Tow = 0.30) 

\ 

\ 

\ 

\ 

\ 

\ 

\ 



Figure 4.— Local skin-friction data (from ref. 2). 



o 

□ 


Sublimation 

Heat transfer (peak heating) 
Prediction (ref. 14) 


10 ' 


R x, tr 


□ 


.013 cm L.E. 
thickness 


0/ 




Jp' 


Zero L.E. thickness 


I o 6 L I I I l I 1 I __ I I I 

3 x I0 6 I0 7 4 x I0 7 


U e / v e , per meter 

Figure 5.— Transition Reynolds numbers for end of transition. 
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Pt.00, N/cm2 

O 394 

□ 666 

O 385 

A 232 

A 452 



22 . 


.20 


p w I8| 
N/cm 2 




B.L. 

survey 

station 


^ ^ i 


^ n 


.14 


12 


.10 


* a ^ ^ ana ^ 


.08 h 


.06 1 


A 


A 


A A 


AA 


A 


.1 .2 


.3 


J_ 

A 


A 


.5 

x / 1 


A 


.6 


A 

. L l4 

.7 


.8 .9 


1.0 


(b) Boundary-layer trips on. 


Figure 6.- Surface-pressure distribution. 
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10 6 10 5 4X10 

P t ( T t )" 7/ f N/cm 2 °K 


Figure 7.— Total temperature probe calibration. 
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T, )e = 1028° K 
H e = 2,590 

O Measured T t 
□ Calculated for Pr = I 



.9 1.0 I 2 


T t /Tt,e 

(a) Profile (a), table 2. 

Figure 8. -Boundary -layer temperature profiles. 



T t) g = 1089° K 
Rq = 5,280 

O Measured Tt 
□ Calculated for Pr = I 



ED 

) 

i I 

.9 1.0 

T t /Tt, e 


(b) Profile (b), table 2 


Figure 8- Continued. 








Re 

T w / T aw 

B.L. trips 

o 

2,590 

.34 

Off 

□ 

5,280 

.32 

On 

o 

4,720 

.45 

Off 

A 

3,960 

.43 

On 

21 

6,910 

.51 

Off 

Q 

8,860 

.51 

On 


a 

T t e « 680° K 

’ /\ 


See table 2 
Linear fairing 


A 


O T t j e « 770° K 


T t e ~ 1050° K 



A 














Subsonic / supersonic 
Fenter , ref. 17 


_U_ 

Ue 


□ 



n 


I 1 


Figure 13.- 


M = 7.4 

Ames 3.5 ft W.T. Wal 



• 2,590 .34 

□ 5,280 .32 

• 4,720 .45 

A 3,960 .43 

A 6,910 .51 

A 8,860 .51 

See table 2 

Solid symbol: B.L. trips off 
Open symbol: B.L. trips on 


J I I I L 

I0 4 

Re 

>file power-law exponent. 






Figure 15.-T-prime reference temperature transformation method (eqs. (5)-(9)); boundary-layer 

trips on. 


omo 




Figure 16- Corresponding-stations transformation method of Coles (eqs. (10)-(1 7)); boundary-layer 


OiDO 



Figure 17.- Corresponding-stations transformation method of Baronti and Libby (eqs. (18)-(24)); 

boundary-layer trips on. 


OcnO 




10 I0 2 I0 3 I0 4 

U T 7 


Figure 18.- Mixing-length transformation method of Van Driest (eqs. (25)-(30)); boundary-layer 

trips on. 



Figure 19. -Mixing-length transformation method of Van Driest; boundary-layer trips off. 
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